Heralded processes on continuous- variable spaces as quantum maps 
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Conditional evolution is crucial for generating non- Gaussian resources for quantum information 
tasks in the continuous variable scenario. However, tools are lacking for a convenient representation 
of heralded processes in terms of quantum maps for continuous variable states, in the same way as 
Wigner functions are able to give a compact description of the quantum state. Here we propose and 
study such a representation, based on the introduction of a suitable transfer function to describe the 
action of a quantum operation on the Wigner function. We also reconstruct the maps of two relevant 
examples of conditional process, that is, noiseless amplification and photon addition, by combining 
experimental data and a detailed physical model. This analysis allows to fully characterize the effect 
of experimental imperfections in their implementations. 



I. INTRODUCTION 

Quantum mechanics is a probabilistic theory. The 
quantum description of any experiment is based on prob- 
ability amplitudes - quantum states and measurements - 
and transformations of such amplitudes - quantum pro- 
cesses. For each of these objects there exist techniques 
to obtain the corresponding mathematical tool from mea- 
sured data: state tomography [L J, , process tomography 
[3l H], and detector tomography [5]. There exist con- 
straints necessary to attribute a near physical meaning to 
abstract mathematics: for instance, a map acting on den- 
sity matrices space corresponding to a physical process 
is normally completely positive (CP). This amounts to 
say that it must send physical states into physical states 
regardless of observing the system by itself or as a part 
of a larger ensemble to which it is de-coupled 0. Most 
of studied maps preserve the norm of the state, but there 
exist notable exceptions: non-trace preserving operations 
arise whenever a measurement on the system is involved. 

An interesting class of such processes involves herald- 
ing: the evolution of a system is considered conditionally 
on the outcome of a measurement on part of the system 
itself [l]. While this is a legitimate operation in classical 
physics, in quantum mechanics non-standard behaviour 
may arise: this is the case, for instance, of anomalous 
weak values [7 . In the context of optical quantum in- 
formation, conditional evolution has found several appli- 
cations for simulating strong nonlinearities at the few- 
photon level. This approach has allowed to build two- 
qubit [3^ and three-qubit quantum logic gates [8^ and 
to generate quantum states with non-Gaussian Wigner 
function [9-13 . Such an evolution is able to induce non- 
Gaussian transformations effective in overcoming existing 
no-go theorems valid for purely Gaussian resources [14]- 
[16]. Successful applications of such processes to commu- 
nication tasks have been demonstrated in several exper- 
iments [T7H2Q]. 

The experimental investigation is relatively at an early 



stage: so far quantum process tomography of non trace- 
preserving maps has been presently implemented only 
in a reduced two-qubits Hilbert space [21^ 22,- Here 
we show, by a detailed physical model, the descrip- 
tion of two conditioned processes which are relevant 
to continuous-variable state manipulation: the noiseless 
amplifier [191 EO] and the single-photon addition JTJ [23] . 
We can derive the expression of the map in the well 
known tensor form, and, as a step further, we illustrate 
a transfer function formalism, which allows to describe 
quantum process directly in the Wigner representation. 
This will stimulate to deepen investigations in this area 
and to develop more sophisticated analytic tools. 

II. QUANTUM MAPS 

Any transformation acting on states needs to satisfy 
some physically-motivated mathematical constraints. In 
the simplest case, a closed system, the evolution of a 
quantum states is described by a unitary operator 

which transforms the input state as p'=Upu\ More 
generally, the system will be able to interact with the 
environment and a representation in terms of a unitary 
won't be sufficient to describe this scenario; however, 
some essential features are retained, in particular the 
output state must be obtained from a linear transfor- 
mation of the input. The proper formalism then adopts 
a generic linear map 8 such that p' — £{p). Similarly 
to the previous expression, this map can be decomposed 
in the incoherent application of a set of Kraus operators 
{Ei} 0: 

£{p) = Y,E,pe\. (1) 

i 

One can note that this expression is similar to the formal- 
ism used for positive-operator valued measure (POVM) 
since a generic transformation can be seen as the appli- 
cation of a unitary operation on a system composed by 
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the input state and the environment, fohowed by a mea- 
surement of the environment for which we do not know 
the outcome. Another expression, more convenient for 
data visuahsation, uses a tensor {Sj^j!^}: 

[^{Q)]l,k = Xl^rr^n,m. (2) 

n,m 

Where the elements of the tensor are given by Ej^j!^ = 

ES\Ei\n){m\El\k) 

These maps can not be completely arbitrary: an essen- 
tial requirement is that they lead physical states in phys- 
ical states. Therefore, these maps have to send positive 
operators into positive operators, and, for deterministic 
processes, require to preserve the trace; so they directly 
give a physical density matrix without need of other op- 
eration. Furthermore, in the majority of the cases, we 
also demand complete positiveness: this amounts to say 
that the evolution must remain physical when the system 
is entangled with a second object. 

A somehow different context arise when considering a 
conditional process : it implies a non-linear evolution of 
the state due to the re-normalisation operation. Indeed, 
these processes often aim to approximate a non-unitary 
linear operator C and transform a pure state \a) into 
^y N{a)C\a) , where N{a) is the normalisation factor, 
which might present a complex dependence on the 
state. And, even if C is actually linear, this linearity 
is shadowed if we co nsider t he physica l inputs, as 
^/N(a)C\a) + ./N{P)C\f3) ^ ^N{a + p)C{\a) + |/3)). 
As a result, in order to represent it as a linear quantum 
map one should use a trace non-preserving map and 
keep the normalisation step for the result of the process. 
Moreover, a map should include one more information 
to describe a conditional process which is the success 
probability. 

In fact, a conditional evolution acts as following: the 
system evolves through a probabilistic device, and we 
only accept those runs when a successful event is flagged 
(Fig. [T]). Clearly, the overall process including both suc- 
cesses and failures can be modelled by a deterministic 
quantum map, and in terms of Kraus decomposition (eq. 
[2| passing from the overall process to the conditional one, 
consist to keep only the subset of the Kraus operators cor- 
responding to the result of the measure used for herald- 
ing. Thus we easily see that the trace of the output state 
gives our additional information: the success probability, 
and a quantum map only need to be trace non increasing 
to correspond to a physical operation [24 . Nevertheless, 
such definition can then be extended to more general pro- 
cesses, as far as the transformation remain linear in the 
quantum state. In particular it can be used with trace 
increasing maps, which, even if they are not physical, are 
very common in theoretical quantum physics and include 
most of the non-unitary operators that are approximated 
by conditional processes. 

The method explained above gives a neat picture of the 
process for discrete-variable systems: for instance, one 



(3) trace-preserving process 




FIG. 1: (a) Trace-preserving quantum operation. The input 
state Q is transformed by the quantum channel 8 in the output 
state S[q]. The trace of the input state is preserved by the 
channel: Tr{^[^]} = 1. (b) Heralded quantum operation. 
The input state q is transformed by the quantum channel 
T in the output state T[q] upon realization of a conditional 
event. The trace of the input state is in general not preserved 
by the channel: Tr{j^[^]} / 1. 

can recognise almost at glance the behaviour of a qubit 
process by inspecting the corresponding tensor. This is 
more complex when dealing with continuous variable sys- 
tems, when often looking at the states as Wigner quasi- 
distribution in the phase space can convey information in 
a more compact and effective way. Therefore, a method 
to represent quantum processes in the Wigner represen- 
tation would be highly desirable. Such an object has been 
proposed for the unitary processes [25l |26] and for Gaus- 
sian operations [15] ; here we give an explicit extension of 
these results to the case of a generic map £. 

For this purpose, we can reason in analogy with the 
probability distribution V{x,p) for physical position and 
momentum of a classical particle. The action of a Marko- 
vian process will modify such distribution via a transfer 
function describing the odds that a parti- 
cle initially in the position will eventually end in 
{x' ^p'). The distribution V'{x' ^p') of the coordinates at 
the end of the process will result from the sum of all these 
elementary displacements: 

V'{x',p') = j dx dpV{x,p)f{x' ,p' ,x,p). (3) 

Hence, we would like to maintain this structure for quan- 
tum processes as well by introducing a suitable transfer 
function by which the input Wigner func- 

tion W{x^p) can be turned into the output W'{x' ^p') by 
the integral transform: 

W'{x',p) = J dxdpW{x,p)fs{x,p\x,p). (4) 

In order to see that this is actually the case, we start 
from the case where only one Kraus operator Ei is 
present; the general result can be obtain by linearity. 
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The Wigner function of the output state then reads: We invoke the completeness relation so to obtain 

W'{x',p') = dve'^^\x' - '^\E^peI\x' + |). (5) 
I 



W'{x',p') = ^ J dpdvdsdte"'P'e'^'-'^PW{^^,p){x' - '^\E,\s){t\El\x' + ^). (6) 
and then, by a variable substitution, the expression for the transfer function associated to the operator Ei: 

I 



which imphes 

J fi{x' ,p' ,x,p)dx'dp' = W^i^{x,p) (8) 
J U{x',p',x,p)dxdp = W^ ^,{x',p') (9) 

Based on the remark that each acts independently, the 
transfer function associated to a generic process reads: 

fs{x,p,x,p) = ^fi{x',p,x,p) (10) 

i 

where each function fi corresponds to a Kraus operator 
with the relation given by Eq. It can be checked 

that using this formula as the definition, and using the 
relation between and Ei^ one arrives to the original 
definition. The transfer function might be a distribution, 
but from Eq. ([7| it appears that it is always real. Also, 
normalisation enforces that 

j fi{x' ,p' ,x,p)dx'dp'dxdp = Tr(£^|£^^) (11) 

= Tt{E,eI) (12) 



Those properties can easily be extended to the whole 
transfer function. In particular, in the case of a deter- 
ministic map we have 



J Mx',p',x,p)dx'dp' = 1, (13) 



whereas for a non-deterministic process the integral of 

W - t {x^p)W{x,p) gives the success probability. 
E- Ei 

We also can express the transfer function in terms of 
the process tensor, starting with the Eq. ([2| rewritten in 
the form 



^(^') = E E ^rrTr {g ■ \n){m\) \l){k\. (14) 

n,m l^k 



We can use the properties of the Wigner functions to 
evaluate the expectation value Tr {g ■ |n)(m|), and derive: 



W'{x',p') = 2n J dxdpW{x,p) (EE^rrW^|n){m|(^'P)^|0(fc|(^>')) ' 



(15) 



\ m,n l^k 



r 



where (^^P) is the Wigner representation of\l){k\. 

Therefore, the quantum operation £ is conveniently rep- 
resented by the transfer function 



fs{x\p\x,p) = 27r^^^^ 



n,m l.k 



{x,p)x 



X Wi 



\l){k\ 



[x'^p'). 



(16) 



The quantum transfer function still bares resemblance 
with Markovian processes. This can be seen by in- 
specting what happens when chaining two processes 
£=£i (8) ^2- Under these circumstances, we obtain the 
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complete transfer function as: 

fs{x' ,p' ,x,p) = J dx" dp" fs^{x' ,p' ,x" ,p")fs^{x" ,p" ,x,p). 

(17) 

which is similar to the Chapman-Kolmogorov equation 
for Markovian processes [26 . This reinforces the view 
that from a classical viewpoint, fs{x' ^p' ^x^p) should be 
interpreted as a a transition probability from {x^p} to 
{x'^p^}. The analogy can not be extended further in the 
quantum domain: in the following, we will illustrate a 
case where fs{x' ^p' ^x^p) can actually take negative val- 
ues. However, we will also show how the temptation 
of establishing a direct quantitative connection between 
nonclassicality and the negative values of fg has to be 
resisted. 



III. DETERMINATION OF THE TRANSFER 
FUNCTION 

In order to determine the expression of the transfer 
function, one could first decompose it in a sum of transfer 
function corresponding to the different heralding events, 
in a similar way as eq. [To] Then, each of those transfer 
function can be constructed by composing, with the use 



of Eq. [T7| of some basic transfer functions corresponding 
to the different elements of the process. The basic trans- 
fer functions can have the same number of input and 
output mode, corresponding to a basic transformation, 
only output modes, corresponding to the introduction of 
an ancilla state, or only input modes, corresponding to a 
measurement. 

The basic transformations can be determined by Eq. [7| 
or by simple considerations. In particular all transforma- 
tion that can be expressed as a coordinate transformation 
have a transfer function composed of Dirac distributions 
which directly come from the coordinate transformation. 
An other interesting basic transformation is the one given 
by a coordinate transformation Mi {i = x^p) for each 
quadrature of the input mode and a quadrature of a vac- 
uum ancilla, followed by a partial trace on the ancilla. If 
the transformation matrix is 

det(Mi) = 1, (19) 

where e = ±1, then the transfer function of this operation 
is: 



f{x',p',x,p) 



exp 



exp 



p - iipp 



(20) 



Table IIIII shows some basics transfer functions obtained 
by those considerations with their tensor form. 

The basic functions with only output or input modes 
are even simpler to determinate. Indeed the first ones 
are exactly the Wigner function of the introduced an- 
cilla, while the second ones are the Wigner function of 
the projector corresponding to the measure multiplied 
by a factor 27r. With those considerations, the transfer 
function can in fact be seen as an extension of the Wigner 
function. 

Finally, the different heralding events are generally the 
ideal heralding and the faulty ones. Consider for instance 
a heralded process when conditioning can be faulty in a 
certain fraction of the total events. We can call Fi the 
correct process, and the failure. The output Qout state 
of the whole process will be a convex combination of 

Qout = i — 5 + (1 - U — 5 ' (21) 

^1 ^2 

where Pi, 2 = Tr[J^i,2(^in)] are the success probability 
used here to normalize the results of both maps. If we 
notice that ^ = Pi /(Pi + P2), we can then infer that the 



transformation which includes both events is given by 

Qout — 5 — [^^) 

where P1+P2 is effectively the trace of ^i(^in)+^2(^m)- 
This amounts to say the correct map is T=T\ + T2 , pro- 
vided that we choose T2 in order to have the correct oc- 
currence probability. This last point is the most crucial, 
and the determination of the good faulty process T2 can 
be complicated. The simplicity of the expression is due 
to the fact that the maps already contains the success 
probability for heralded process. Nevertheless we should 
notice that if the false heralding comes from a noisy mode 
in the input one should add (at least) a supplementary 
input mode to the map. 



IV. EXAMPLE 1: THE NOISELESS AMPLIFIER 

We now inspect two important quantum processes with 
the formalism of quantum maps: we will be able to high- 
light clear signatures of nonclassicality, and observe how 
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TABLE I: Tensor process and transfer functions of some basic transformations 



Transformation 


tensor 


transfer function 


Identity 




6{x — x')5{p — p) 


Phase rotation 




5{x — cos Ox' + sin^p^)^(p — cos Op' — sin^x^) 


Displacement 


Vm!n!/!A;!ei i^^^o Z^j=o J j 

(l-i)\(k-j)\ ^ ^ 


o\x — X ^ y 1)xc\ol))0\p — p ^ y IVvc^yoL)) 


Squeezing 


(m!n!o3/2 cosn(^rj 2^ .^^ i^.^Q A 


8(x — e^x)6(p — e~^p) 


Beam splitter 


/mi!m2!ni!n2! Y^^l V^^l / 1^1 + ^1— 

x('^)( ) 

\% J \m\—%J \ J / yni—j/ 
^^2i+2j + l2 + k2—mx—n-i^lx+kx+m^+n-i—2i — 2j 

V n 1 1 \ 1 n 1 / 1 / 
^'-'mi+m2,ti+t2 ^n\-\-n2 ,k\-\-k2 


6{xi - tx'i + rx'2)5(j)i - tp'i + rps) 
x5{x2 - tx'2 - rx'i)5{p2 - tp'2 - rp'i) 


down- conversion 


/ mi !m2 'tt-i !n2 ' (P"!) 

A/ /i I/,-. I t-i I Z^^ I ('r7^1 +71,1 +Zo +fco ") /2 + 1 
y ti 112 •"'I •"'2 • (jf V '"'1 "'1 ^^"-2 / ''T^-^ 

^ ^n2 —ni,k2 — ki ^m2 —'m-i,l2—l\ 
Y sr^li (n2+^)!(m2+j)! 
^ Z^i=0 Z^j=0 (ni-fei+'i)!(mi-Zi+j)! 
(_l)*+igfci+^i-i-j((^_l)i+j-(fci+ii)/2 

i\(kl-i)\j\(h-j)\ 


^\\/y^^ ^ yy ^^2 ^1 ypi yy ^P2 pij 
x5{y/gx2 + \/9- Ixi - X2)6{y^p2 - Vg - Ipi -P2) 


Attenuation 






V i!fc! (m-0! Om-l,n-k 


7r(l-r7) ^""J^ I 1-7] 1-7] ) 


Parametric amplification 


/l\k\2^ + m 2^-"^ gk + (m-n)/2 


_L_exp(-(£:^-(i^:^) 


V n\m\2^ + '^ (fc-m)! ((^ + i)fc + n + l Oi-n,fc-m 





they degrade under experimental conditions. As an in- 
teresting feature, we will be able to capture such non- 
classical aspects at a glance. Our analysis first concerns 
the noiseless amplifier {C = where h is the number 
operator and ^ > 1) [l9l|20]. We do not adopt a "black 
box" approach, rather a model of the process is used so 
to arrive to a description in term of generalised maps, 
also in the case when all the imperfections are taken into 
account. 

Our device (Fig. |2| is the teleportation-based ampli- 
fier proposed in Ref. [27 : its working principle is to 
use a non-maximally entangled resource - a single pho- 
ton split on an asymmetric beam splitter (A-BS) - to 
perform the teleportation of a coherent state \a). The 
analogue of the Bell-state measurement consists of su- 
perposing the reflected portion of the single photon with 
the input state on a symmetric beam-splitter (S-BS), and 
perform photon counting at the outputs. Successful runs 
are heralded by the presence of a single photon on one 
output and the vacuum on the other. This operation 
produces an output state in the form N{a) (|0) + ga\l))^ 
where ^ is a gain factor dete rmined by the reflexion R of 
the A-BS, g = y^(1 — R)/R. For weak input intensities 
II a II < 0.1, this truncated expansion is a good approxi- 
mation of the amplified state \ga). Figjs] (a) shows the 
elements J~'^jJ^ of the corresponding map, which would 
normally describe the population transfer among Fock 
states. In this case, instead, we notice the enhancement 
of the single photon component and the suppression of 
all the higher-order terms. The complete process is a 
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FIG. 2: Layout of the implementation of the noiseless am- 
plifier [19]. A single photon is conditionally generated upon 
detection of a single photon on detector Do- After splitting 
in an asymmetric beam-splitter (A-BS), the single photon 
is mixed with the input coherent state \a) in a symmetric 
beam-splitter (S-BS). The noiseless amplification process oc- 
curs conditionally to the detection of a photon on detector 
Di. The beam-splitter operations are performed exploiting 
polarization (double sided arrows in the figure). Inset: fideli- 
ties between the experimental density matrices 19 and the 
prediction of the model. 



truncated form C = ^ 0(n), where 0(n) is 1 for n < 1 
and zero otherwise. 

Several departures from the ideal behaviour prevent 
from matching these simple predictions in the experi- 
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(a) 



(b) 





(c) 



(d) 



FIG. 3: (a) Diagonal elements J^^lJ^ of the ideal truncated 



noiseless amplifier process, (b) Diagonal elements J^^jJ^ with 
non-unit detection efficiency of the APD Di (/i = 0.11) and 
lack of photon-number resolution, (c) Diagonal elements 
-^k^k^ with non-ideal generation of the single-photon state 



(S = 1.089) (d) Diagonal elements J^^ 
perimental imperfections. 



■m,m 
k 



including both ex- 



ment, and a more refined description is then necessary. 
One of the main limitations is represented by single- 
photon detection. While the apparatus is quite robust 
against limited efficiency [19j |27], it is nevertheless af- 
fected by the lack of photon-number resolution. The 
APD Di, which heralds the successful events of the am- 
plification process, will give a click each time that some 
light is absorbed, irrespectively on the energy. This 
causes triggering events which do not originate from sin- 
gle photon on Di , which result in a transfer of population 
from higher-energy states to the one-photon Fock state, 
as it appears in Fig. [3](b), where the corresponding ten- 
sor is shown. 

Single photons are produced by down-conversion in a 
non-linear crystal: whenever an avalanche photo-diode 
(APD) Do detects the presence of one photon, it heralds 
the twin photon on the correlated mode. The probabilis- 
tic nature of the emission allows for multiple-pair gener- 
ation, which the APD is not able to discriminate from 
single-pair events. The output state will not be a pure 
single photon, but will present contributions from higher 
order terms. Furthermore, one needs to consider that the 
matching of the pump field with the observed modes will 
not be perfect. This results in excess noise in both the 
conditioning and signal modes, spoiling even more the 
quality of our single photon state. This can be assessed 
in the experiment by measuring the quadrature distri- 
butions , and is described by a parameter S [28] |29] , 
ranging from (5 = 2 for a pure heralded single-photon 



state to ^ = for a thermal state. The effect on the pro- 
cess is illustrated in Figjs] (c) showing how higher-order 
number from multiple-pair emission can end up being 
populated. 

Considering both imperfection sources provides an ex- 
haustive model of our experiment: its accuracy can be 
checked by calculating the fidelities between the exper- 
imental density matrices [19 and the prediction of the 
model (right inset of Figj2|, showing an average figure of 
~ 99.5%, for weak intensities ||a|| < 1. The results are 
summarised in Fig. |3|d): the two mechanisms take place 
independently and cause a lesser gain than expected, and 
the presence of noise in the amplified states. 

The transfer function can be decomposed in a correctly 
heralded and a faulty one. For the correctly heralded 
transfer function, one could start with the beamsplitters: 
the asymmetric one is obtained by composing, with the 



use of Eq. 17 , the beamsplitter transfer function with the 
Wigner function of the experimental single photon and 
the vacuum. On the same way, the transfer function for 
the part containing the symmetric beamsplitter is deter- 
mined by composing the beamsplitter transfer function 
with the transfer function of the APD (Wigner function 
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FIG. 4: Contour plots of fj^{r\r,0) for the noiseless amplifi- 
cation process as a function of {r , r) for different values of 0. 
(a) Low APD detection efficiency /i and ideal generation of 
the single photon (^ = 2). (b) Low APD detection efficiency 
/i and non-ideal generation of the single photon (S = 1.089). 
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of the projection operator multiplied by 27r composed 
with an attenuation transfer function) on one output and 
an attenuation (modelising the mode-matching on one 
input), and then tracing on the remaining output mode. 
The final correctly heralded transfer function is obtained 
by composing those two transfer function. The faulty 
transfer function is simply determined by composing an 
attenuation transfer function (taking into account the 
mode-matching and the symmetric beamsplitter) with 
the APD transfer function. 

While an inspection in the Fock basis can be informa- 
tive, it does not lead to the most natural description of a 
process for continuous- variable states; also, from a practi- 
cal point, it might be cumbersome to verify some proper- 
ties such as the Gaussianity of the process, or its nonclas- 
sicality from the expression of the T tensor. On this pur- 
pose, a useful approach consists in inspecting the trend 
of the associated transfer function: while this object gen- 
erally acts on pairs of two-dimensional vectors r — (x^p) 
and = {x' ^p')^ for phase-invariant processes - as it is 
the case for the noiseless amplifier - the transfer function 
can only depend on r = Y^^x^^Tp^, r' = x'^ ^ p"^-, and 
6 = cos~^ ( ). The transfer function for the noise- 

less amplifier is presented in Figj4j comparing the cases 
when an ideal single photon is used as ancilla and with 
the actual resources. Non-classical features are clear in 
the ideal limit, in which negative values appears around 
r'=0. However, in experimental conditions, these sig- 
natures are smoothed by the imperfections of the set-up, 
though there remains negative region. In more details, we 
observe how low values of the transfer function for high 
correspond to the saturation of the amplifier, i.e. the 
impossibility of having more than one photon at the out- 
put. The negative peak determines the non- Gaussianity 
of the output states by causing a small negative region 
in the Wigner function of the output state; nevertheless, 
this feature vanishes rather quickly with the imperfection 
and is not visible in realistic output state, as it appears 
in FigJH We can notice that the region for small r and r' 
is quite different around 6 = and ^ = tt, this is due to 
the fact that the amplifier keeps the phase of the " small" 
states, whereas for bigger value of r and r' the trans- 
fer function is almost independent of since the higher 
photon-number terms in the input state often trigger the 
heralding leading to a single photon (with losses) in the 
output state. A last remark concerns the increasing peak 
with r and the different scales between the two maps orig- 
inating from the variation of the success probability. 



V. EXAMPLE 2: PHOTON ADDITION 

The extreme negative value of a Wigner function can 
be used as a quantifier of its non-classicality [23]. How- 
ever, an intuitive extension of such a reasoning to the 
transfer functions would be severely affected by the prob- 
abilistic character of the process itself. Here we illustrate 




FIG. 5: Layout of the single-photon addition experiment [23] . 
A single photon is conditionally added upon detection of a 
single photon on detector Do- As above, the beam-splitters 
exploit polarization (double sided arrows in the figure). Inset: 
fidelities between the experimental density matrices [23 and 
the prediction of the model. 



these considerations in a second example: the single- 
photon addition (C=d^). As above, our description is 
mediated by a model of the physical process. 

In our implementation, photon addition is achieved by 
feeding the input state in an optical parametric amplifier 
(OPA) driven at low gain g=cosh.^ Xi where x measures 
the non-linear interaction strength and it is proportional 
to the pump intensity. To the first order, this process 
adds a photon pair shared by the signal mode, and a 
correlated mode, on which an APD I^o is placed; due 
to the nondeterministc nature of the process, the suc- 
cessful events are triggered by a detection event from 
Do (Fig. [5|. This method has been introduced in [TT], 
and then adopted for tests of the commutation rules 
j3Q| [3T] , and the analysis of non-classicality [32^ and non- 
Gaussianity [23] . 

The main source of noise here can be identified in the 
imperfect matching between the pump and the signal 
modes, which results in a parasite gain h= cosh^ jx- An 
estimate for these two gains, taken from a fit of their non- 
Gaussianity [23 , is x=0.105, and 7=0.425 [23,. A third 
imperfection arises from the fact that spurious events 
might happen at I^o, due either to dark counts or clicks 
originating from non- matching modes. The average fi- 
delity between modelled and the reconstructed states us- 
ing coherent states as inputs is satisfactory, although the 
data might be affected by some extra noise likely due 
to low-frequency fiuctuations of the average level of the 
homodyne current. 

As it appears from the comparison of Figs. and b, 
the gain x is chosen to be sufficiently low so that two-pair 
events are not significant: the transfer of population by 
more than one photon is low. On the other hand, the 
effect of the parasite gain seems as important: the sheer 
effect is the presence of uncorrelated clicks at I^o that 
leave the state unchanged. This corresponds to the diag- 
onal terms in Fig.jGj), considered in the limit of extremely 
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low gain x ^ 0. The overall process simply results in the 
presence of these two imperfections. For the sake of sim- 
plicity we have only considered the case of low detection 
efficiency at I^o (Fig.|6]i). In this example, the adoption 
of the quantum map formalism reveals to be particularly 
clear and useful for the analysis of the process: not only 
it conffims our intuition about the behaviour of para- 
site processes, but also give us a way of quantifying their 
effect in a way that does not depend on the particular 
input. 

The correctly heralded transfer function is easy to de- 
termine: it is only a composition of the parametric down- 
conversion transfer function with the Wigner function of 
vacuum on one input and the APD transfer function and 
attenuation transfer function on the outputs. As for the 
preceding example, we can use the radial symmetry of the 
process to express the transfer function in the simplest 
form fjr{r'^r^9). In Fig j?] we show the transfer function 
,p as a function of (r',r), as this contains most 
information about the physics: there we compare the case 
of low gain r and photon-number discrimination, and the 
full model of our experiment. For both cases, the trans- 
fer function has non-zero values only around r' = r and 
for = 0: this indicates that the amplitude and phase 
are mostly unchanged by the process. The increase of 
the positive peak with the amplitude corresponds to the 
growth of the success rate with the number of photons, 
and is more visible in the second graph because of the in- 
ability to discriminate the photon number. Finally, the 




(a) (b) 



FIG. 7: Plots of fj^{r,r,0) for the photon addition process 
as a function of (r\ r) for ^ = 0. (a) With ideal OPA (driven 
at x=0-105) and photon counter, (b) With a parasitic gain 
7=0.425 and APD. 




-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 
r'-r r'-r 



(a) (b) 

FIG. 8: Plots of fj^{r,r,0) for the photon addition process 
as a function of r' — r. The plain line correspond to the 
experimental conditions ( 7=0.425 and APD), and the dashed 
line to the ideal case, (a) for ^ = and r + = 2 (b) for 
6* = and r + = 20 





negative peak introduce a negative part in the resulting 
Wigner function and is a sign of the non-gaussianity of 
the process. We can note that the difference between the 
two scales (Fig. [8| are only due to the differences in the 
success rate, even for the negative peak. It implies that 
the size of the negative peak can not be readily used to 
quantify the quant umness of a map. 



(a) 



(b) 



VI. CONCLUSION 





(d) 



FIG. 6: (a) Diagonal elements J~'^lJ^ of the ideal photon ad- 
dition process, (b) Diagonal elements J^^lJ^ for the case of a 
conditioned OPA driven at x=0.105. (c) Diagonal elements 
^k^k^ with a parasitic gain 7=0.425 and very low gain (d) 
Diagonal elements J-'JTlJ^ including both experimental imper- 
fections. 



We have inspected two important processes for 
continuous-variable states with the formalism of quan- 
tum maps: these can convey interesting physical informa- 
tion about the process independently on the state. The 
adoption of a description in terms of transfer functions 
offers a compact, insightful view of the process, along the 
same lines of what happens with the Wigner function for 
quantum states. We have applied this method to the de- 
scription of the realistic operation of a noiseless amplifier, 
and of a photon-adder, evidenciating how experimental 
imperfections shape the features of the transfer function. 
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EU project ANR ERA-Net CHISTERA HIPERCOM. 
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